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Stimulated by the pioneering work of Peres |l],|2| on state sensitivity to small changes 
in the Hamiltonian, several new developments have taken the study of quantum state sen- 
sitivity to the experimental realm. Schack and Caves 0] have analyzed the sensitivity to 
perturbations as the amount of information about a perturbing environment that is needed 
to keep entropy from increasing. Their approach links the study of quantum sensitivity to 
statistical mechanics and to quantum information theory. Recently, Gardiner et al. Q have 
proposed a detailed experimental set-up of an ion trap to study quantum state sensitivity to 
changes in the Hamiltonian. A different approach has been taken by Ballentine and Zibin || 
that study the emergence of classical state sensitivity from quantum theory in computational 
time reversal. The difficulty of directly tackling the problem of quantum state sensitivity to 
initial conditions stems from the following reason. Classical sensitivity to initial conditions 
is made precise when measured by Lyapunov characteristic exponents || 



with d c (x t , y t ) the Euclidean distance between two phase space vectors x t and yt. There 
is classical sensitivity to initial conditions when A c (xo) > and the trajectory starting at 
x is said to be unstable. A direct attempt to define the quantum analog of the classical 
Lyapunov characteristic exponent in (fj) substitutes the phase space distance between two 
classical trajectories by the Hilbert-space TC distance between two close quantum vectors 
giving 



a zero quantum Lyapunov exponent for all states as any two states do not separate at all with 



d H mt))Av(t))) = \\\m)-W))\\ = IIIV'(O)) -kp(0))|| and ll-M = /<T>- According to 



conditions can be found. 

However, there are three strong objections against the use of as a measure of state 
sensitivity to initial conditions and instability. Firstly, the quantum mechanical state 
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(0) all possible quantum states are stable and no sensitivity of quantum states to initial 
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space is the complex projective Hilbert space CV n that has a curvature not present in 
the complex Hilbert space 7iJ +1 . It is a nontrivial Kahler manifold with a symplectic 
form and an associated Riemannian metric ||. Secondly, equation (|2|) is obtained by 
substituting the classical phase space trajectory distance in ([I]) by a quantum distance. 
This is not a mere technical error but a conceptual one as the closest classical object 
to a quantum state is a Liouville density, not a single trajectory HfJJ^]. Recall for ex- 
ample that in the classical limit the dynamical equation for the Wigner function or for 
the coherent state representation of the quantum density reduce to the Liouville equa- 
tion 0. A definition analogous to (fj) for Liouville densities gives Xc [p(0)] = with 
dc (pi(t),p 2 (t)) = \J J dqdpp\ /2 (q, p, t)pl /2 (q, p, t) = d c (pi(0), p 2 (0)) by Koopman's theorem. 
The fact that Xc(p(t)) — for all possible Liouville densities tells us that we cannot know 
from (|2|) if quantum states have sensitivity as this type of measure for classical Liouville 
densities is insensitive to classical-mechanical instability. Thirdly, unlike the classical tra- 
jectory distance, d n is a bounded metric, thus precluding the exponential divergence. In this 
paper we propose to overcome the above three objections in three steps to provide a working 
scheme to classify quantum and classical-statistical states by their sensitivity and instability 
in projective space V n . This sensitivity and instability are shown to be measurable by the 
different transition probabilities of two nearby states to a third fixed state. 

Step 1. A classical density p(x) can be mapped into a real Hilbert space rC£ l by taking its 
square root, ip c {x) = \J p(x). The expectation of a classical operator F is (F) = (Fipf) / (-0^) , 
which means that the physical space is not 7i^ +1 but the space of equivalence classes obtained 
by identifying ip c ~ \ip c for any A G R — {0}. This physical space is the real projective 
space RP n . Similarly, in quantum mechanics the true physical space is not the complex 
Hilbert space Ti.Q +1 but the space of equivalent classes obtained by identification of vectors 
ipQ ~ XipQ for any A G C — {0} , the complex projective space CV n . Classical and quantum 
states in V n will be written as ip. The distance between two points in V n (real or complex) 
can be defined as 
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(3) 



which is the length of the geodesic connecting the two points in V n as measured by the 
Fubini-Study metric ||[| . 

Step2. Now we are faced with the problem of finding relevant distances to characterize 
the motion of classical-statistical and quantum systems by their sensitivity and stability 
properties. According to the quantum expectation rule f2|, the probability that a quantum 
state prepared in state ips will pass succesfully the test for the state ipR is (ips V'R 



This 

shows that the distances of two arbitrarily close states ips and ips + Sip to a third fixed state 



ipn are experimentally relevant as from ([3]) (ips ipR 



for 



cos 2 (d-p ips,ipn /2^ and similarly 
are different and 



ips + 5(pj . These distances dp ips(t),ipR and d-p ips(t) + S(p(t),ipR 
time-dependent and their difference is expected to behave differently in stable and unstable 
systems. 

Step 3. The distance between two states in (|3|) is bounded, d-p G [0, ir] , which obviously 
precludes exponential divergence. To understand how to study unstable systems in terms of 
a bounded metric we express the unbounded phase space Euclidean distance, d c G [0, oo) , 
in terms of a bounded distance, say d\ G [0, it] , that preserves the topology (d c is iff d c is 
0) as |0| 



d\ (x, y) = ixd c (x, y) / (1 + d c (x, y)) 



(4) 



The classical-mechanical Lyapunov exponent in (|TJ) can then be expressed in terms of the 
classical bounded metric in ffl) as 



A c (x ) = lim lim - In (A c (x t , y t ) /A c (x , y )) 



(5) 



with A c (x, y) the divergence function that is the unbounded phase space distance expressed 
in terms of the bounded one as 



A c (x,y)EE^ (x,y)/(7r-d$ (x,y)). 



(6) 
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We are now ready to overcome the three objections to the definition in (0) and propose an 
analogue of the classical-mechanical Lyapunov exponent in V n (that we name P-Lyapunov 
exponent). Let ips(t) and ipsif) + 8(p(t) be two arbitrarily close system states and ipR a 
reference state on the path ips(t), ipR — V , 5( r ) • The P-Lyapunov exponent is obtained by 
the different distances of the two arbitrarily close system states to the reference point of the 
form 
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with A 



ips{t) the divergence function for the distance between the system state ips{t) and 



the reference state if) a 



A 



d 



V 



M^s{t) /(ir -dp ip R ,ip s (t) 



(8) 



and dp the bounded distance between states in (|3]). There is sensitivity to the initial con- 
dition ips(fy of degree A^ when Ap° > and the path ipsit) is sa id to be unstable. The 
calculation of Ap can be obtained with a single path by choosing ipsif) +6<p(t) = tps(t + 5t). 

We have obtained the definition of P-Lyapunov exponent in (0) by overcoming the three 
objections to the definition in fl2|). We now put this definition to the test. The approach 
taken to find the 'P-Lyapunov exponent has been kinematic and it is expected to be valid 
for classical-statistical and quantum states. First we show its validity for classical-statistical 
states and its equivalence to the classical-mechanical Lyapunov exponent for several exam- 
ples. Consider the classical unstable map <&(x) = rx with r > 1, for which the classical- 
mechanical Lyapunov characteristic exponent is A c (xo) = lnr. After n steps the classical 
classical density evolves to p n (x) = P n p(x) = p(x/r n ) jr n . Take as initial system state a 
square- looking function -0s (x, 0) - 



' Po(x) = y Q(x) — Q(x — b) and the Heaviside func- 
tion . The distance between the system state and the reference state x/jr = ips{fy is found 
to be dp = 2 arccos |r~ n//2 |. Its corresponding divergence function ([8]) ultimately diverges 
exponentially and the V— Lyapunov exponent in (0) can be calculated to be 
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ln(r)/2 = A c (x )/2, 
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(9) 



half the classical- mechanical Lyapunov exponent. There is then sensitivity to initial condi- 
tions for classical statistical states and they are classified as unstable in the same manner 
than using classical- mechanical Lyapunov exponents for all r. Figure 1(a) shows the con- 
vergence of \ l v to for different values of r. 

A more sophisticated map is the r-adic transformation = rx (modi) with r > 1, 
which is an exact transformation || (therefore mixing) with classical-mechanical Lyapunov 
exponent A c (xo) = lnr. A necessary feature added to instability for mixing systems is bound- 
edness that implies, as in the classical-mechanical trajectory case, that the distance between 
the system state and the reference state is saturated at the value 2 arccos | ((ips) (^r)) I , but 
before this saturation takes place the P-Lyapunov exponent can be extracted to be arbitrar- 
ily close to ln(r)/2 = A c (x )/2 by choosing a sufficienty localized initial classical-statistical 
state. A classical transformation related to the r-adic map of interest for us in the later 
discussion of quantum states is the baker's map, a map of the unit square onto itself of the 
form 

*{x,y) = (2x-[2x],(y + [2y])/2) (10) 

where [•] indicates integer part. This map is area-preserving and mixing with classical- 
mechanical Lyapunov exponent In 2. Following the same steps as in the r-adic map, we find 
the same result with r = 2. 

We have obtained classical-statistical state sensitivity to initial conditions and insta- 
bility for classical-mechanical unstable and mixing maps. Encouraged by these results we 
now proceed to study quantum systems. To test our definition of V— Lyapunov exponent for 
quantum states we first choose simple analytic cases, the parabolic barrier and the harmonic 
oscillator. All states of the harmonic oscillator (see Figure 1(b) for a Gaussian function) 
are stable, Ap = 0. For the parabolic barrier, with potential V(x) = —moj 2 x 2 /2, and 
classical-mechanical Lyapunov exponent A c (xo) = u there are stable and unstable quan- 
tum states. The most stable ones are the eigenstates cj) a and in fact any discrete sum 
of eigenstates ips = J2n c n<fin is stable. On the other hand, states that use the contin- 
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uous spectrum are unstable. For example, the distance between a Gaussian initial state 
on top of the barrier, ips{ x ) — ex P {—mujQX 2 /2h) , and the reference state is of the form 
2 arccos | (cosh 2 ut + (u/ujq — ujq/uj) 2 sinh 2 out) ^ |. Its corresponding divergence function 
(§) ultimately diverges exponentially with time and the "P-Lyapunov exponent in (0) is 
found to be 



\ oo 
A-p 



MO) = cu/2 = A c (x )/2. (11) 



Therefore, in contrast to the starting definition (0), there is quantum state sensitivity to 
initial conditions and quantum unstable paths. Figure 1(b) shows the convergence of \ l v to 
the "P-Lyapunov exponent Xp for two states with different values of u. 

We now consider the added complication of classically mixing systems for which we 
will use the classical baker map in ([U]). Quantization of a classical map is not unique 
and we consider two quantum baker maps to illustrate two different types of behaviour. 
First we understand the map as defined in configuration space, as did Mendes [[11]] for 
the Arnol'd cat map, for which a quantum theory can be defined by the unitary Floquet 
operator U such that ips{x,y) — > %l)' s (x,y) = ips ($ _1 (^, y)) = Uips( x ,y) with $ in ( [TDD 
and p Xt y = —ihV Xt y the momentum. For this quantum map the situation is analogous to 
that of classical-statistical states under the baker map and one then finds state sensitivity 
and unstable quantum paths in the same way. The second quantization procedure that we 
consider understands the classical baker map in position-momentum space and we adopt the 
Balazs-Voros-Saraceno [0)0 quantization procedure. The new characteristic that interests 



us is the discreteness of the spectrum. As we have already found for the parabolic barrier, it 
is possible to have stable states made of a discrete sum of eigenstates |15|,^4|,|16|,|17| even when 



localized onto classically unstable structures [{18 ]. However, this does not pose a problem for 



the correspondence principle as for systems that appear classical to us the recurrence times 



of the quantum states are longer than the age of the Universe ||19|| . In typical numerical 
and experimental results the recurrence times due to the discreteness of the spectrum might 
be longer than that due to classical boundedness by the energy hypersurface. We relax 



the mathematical convenience of considering asymptotic definitions and consider coverged 
values for experimentally relevant times. Observe in Figure 1 that nonasymptotic values 
already give a classification of states by their sensitivity and instability for finite times. We 
now proceed to study the Balazs-Voros-Saraceno quantum map. Position and momentum 
are discretized qj = Pj = (j + 1/2) /N for j = 0, iV — 1 with iV even. The transformation 
between the position basis {1^)} and the momentum basis is given by (Gn)^ = (Pi\ Qj) 

and the baker map is then defined as 



B = G N l 



f G N/2 ^ 
i Gn/2 j 



(12) 



where the matrices Gn/2 map position kets on the lower (upper) half of the position interval 
to the lower (upper) part of the momentum interval and G N X maps back to the position rep- 
resentation. We choose as initial quantum state ips (Qj) = exp (q Q — qj) 2 /2a + ip qj/a^j 
with a = 10 4 and qo = Po = 0.003. Figure 2(a) shows the finite-time V -Lyapunov exponent 
for iV = 1800 and times shorter than the bounded time i& ~ 8 and it is found to be close to 
half the classical-mechanical Lyapunov exponent ln(2)/2 as in the previous examples. 

The notion of V -Lyapunov exponent introduced in this paper can be measured in many 
types of experimental set-ups. It can be measured in optical filtering experiments by the 
different reduction of intensity of two nearby states. In molecular experiments the distance 
between the system state and the reference state can be obtained from the Fourier transform 
of the Franck-Condon spectrum. Very recently, Gardiner et al. @] have proposed an ion trap 
set-up as a practical realization of a delta-kicked harmonic oscillator to measure the Peres || 
overlap 0(t) = (ips\ U\(t)U 2 (t) \ips) with U\ and U 2 two close unitary evolution operators 
to study the state sensitivity to changes in the Hamiltonian. We can reinterpret their results 
as sensitivity to initial conditions in the following manner. Consider the two nearby vectors 
IV's(O) = Ui (t)U 2 {t) (^(O)) and \ij)s{t + dt)) . Their distances from the reference vector 
\iPr) — |^s(0)) are given in terms of the Peres overlap by 2arccosO(t). Figure 2(b) shows 
the divergence function for the two states of Fig.4(b) and (d) of [§J, respectively. Figure 2(b) 
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also shows the V -Lyapunov exponent and we find that one path is unstable with Xp rs 0.017 
and the other stable with ~ 0. 

The classification performed by the quantum P-Lyapunov exponent will be in general 
different from the classical-mechanical one, but they must coincide in the classical limit. To 
see this first note that the Wigner representation or coherent state representation of the 
quantum state reduce to the classical Liouville density in the classical limit 0. Moreover, 
the classification performed by the classical-statistical (Liouville) P-Lyapunov exponent co- 
incides with the one given by classical-mechanical Lyapunov exponents. To see this note 
that the distance between two classical densities given by (|3|) when the underlying trajec- 
tories obey x(i) = x(0) exp(A c t), with A c the classical-mechanical Lyapunov exponent, is 
obtained to be dp ips^R — arccos (exp (— A c t)) and the classical Liouville P-Lyapunov 
exponent reduces to X-p ^s(O) = A c /2 . This property opens the possibility of studying in 
detail the classicalization of instability, most interestingly by the calculation of the quan- 
tum "P-Lyapunov exponent when adding classical or quantum noise and with increasing 
dimensionality of the system. 

To conclude, a successful measure of quantum and classical-statistical state sensitivity 
to initial conditions and instability has been obtained. The classification given for quan- 
tum states reduces in the classical limit to that performed by classical-mechanical Lyapunov 
exponents. This classification has been studied in several examples, including recently pro- 
posed ion trap states. 

J.I. Cirac, M.S. Child, S. Gardiner, R. Schack, E. Sjoqvist and P. Zoller are acknowledged 
for useful discussions. A Marie Curie fellowship is also acknowledged. 



9 



REFERENCES 



[1] A. Peres, Phys. Rev. A 30, 1610 (1984) 

[2] A. Peres, Quantum Theory: Concepts and Methods (Kluwer Academic, Dordrecht, 1993) 

[3] R. Schack and CM. Caves, Phys. Rev. Lett. 69, 3416 (1992); 71, 525 (1993); Phys. 
Rev. E 53, 3387 (1996);Phys. Rev. E 53, 3257 (1996). 

[4] S.A. Gardiner, J.I. Cirac and P. Zoller, Phys. Rev. Lett. 79 (1997) 

[5] L.E. Ballentine and J.P Zibin, Phys. Rev. A 54, 3813 (1996) 

[6] A. Lasota and M.C. Mackey, Chaos, Fractals and Noise, Stochastic Aspects of Dynamics 
( Springer- Verlag ,1994) 

[7] We use the word distance for the classical-mechanical case although the phase space 
has no metric structure but a symplectic structure. However, the metric can be fixed by 
considering the classical limit of the quantum metric for minimum uncertainty states as 
in J. Anandan , Phys. Lett. A 147, 3 (1990). In any case, for an arbitrary fixing of the 
metric an infinitesimal distance between trajectories will ultimately grow exponentially 
for chaotic systems. 

[8] J.E. Marsden and T.S. Ratiu Introduction to Mechanics and Symmetry (Springer- 
Verlag, 1994) 

[9] J. Anandan and Y. Aharonov, Phys. Rev. Lett. 65, 1697 (1990) 

[10] J.R. Giles, Introduction to the Analysis of Metric Spaces (Cambridge University Press, 
1987) 

[11] R.V. Mendes, Phys. Lett. A 171, 253 (1992) 

[12] N.L. Balazs and A. Voros, Ann. Phys. (N.Y.) 190, 1 (1989) 

[13] M. Saraceno, Ann. Phys. (N.Y.) 199, 37 (1990) 



10 



[14] T. Hogg and B.A. Huberman, Phys. Rev. Lett. 48, 711 (1982) 

[15] P. Bocchieri and A. Loinger, Phys. Rev. 107, 337 (1957) 

[16] R.C. Brown and R.E. Wyatt, Phys. Rev. Lett 57, 1 (1986) 

[17] G. Radons and R.E. Prange, Phys. Rev. Lett. 61, 1691 (1988) 

[18] G. G. de Polavieja, F. Borondo and R. Benito, Phys. Rev. Lett. 73, 1613 (1994) 

[19] A. Peres, Phys. Rev, Lett. 49, 1118 (1982) 



11 



FIGURE CAPTIONS 



Figure 1. Classification of states by their sensitivity to initial conditions as measured 
by the V -Lyapunov exponent in (0) for (a) classical-statistical states under the unstable 
map = rx with r = 2, 3 and 5 converging to half the classical-mechanical Lyapunov 
exponent ln(r)/2 and (b) for quantum states in the harmonic oscillator and the parabolic 
barrier with uj = 2 and 5 converging to half the classical Lyapunov exponent 0, 1 and 2.5, 
respectively. 

Figure 2. Classification of quantum states with discrete spectrum by their finite-time 
V -Lyapunov exponent. Pannel (a) shows the finite-time V -Lyapunov exponent for the 
quantum baker map close to ln(2)/2. The probability density is shown as an inset. Pannel 
(b) shows the divergence function (H) for a state centered on a classically regular region and a 
state centered on a classically chaotic region of a delta-kicked harmonic oscillator, indicated 
by (1) and (2) respectively. The inset shows their corresponding finite-time V -Lyapunov 
exponent that converges approximately to and 0.017. The data for (b) has been obtained 
from the experimental proposal of Gardiner et al. j|] on an ion trap set-up. 
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